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tH- , Abstract 

CN 

[~^ ' In order to benchmark the 3-diniensional calculation of the atmospheric neutrino 

^D ■ flux based on the FLUKA Monte Carlo code, muon fluxes in the atmosphere have 

been computed and compared with data taken by the CAPRICE94 experiment at 

ground level and at different altitudes in the atmosphere. For this purpose only 

two additions have been introduced with respect to the neutrino flux calculation: 

, the specific solar modulation corresponding to the period of data taking and the 

Q I bending of charged particles in the atmosphere. Results are in good agreement with 

r^ • experimental data, although improvements in the model are possible. At this level, 

however, it is not possible to disentangle the interplay between the primary flux and 

the interaction model. 
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1 Introduction 

Recently, a new atmospheric neutrino flux calculation, based on the FLUKA Monte Carlo 
model [0], has been proposed 0. This new calculation has been produced with the aim to 
introduce a high level of refinement in the modeling of all aspects of shower calculations in 
atmosphere, from particle production to particle transport. The need for a 3-Dimensional 
approach has also been demonstrated. For example, considering as reference the Bartol 
flux calculation p], and accounting for the differences in the 1-D and 3-D approaches, the 
FLUKA results are, for E^ > 1 GeV, ~15% lower. This has been obtained for the same 
input primary spectrum and the reduction is present also at energies far from effects due 
to geomagnetic field. This was attributed to the differences between particle production 
models|^. An almost identical difference is present with respect to the Honda et al. 
simulationPI (HKKM in the following), but in this case the input primary spectrum is 
not the same and is believed to be the main source of the discrepancy. 
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Since at present there is not yet a firm experimental determination of tlie absolute 
neutrino flux, one of the possible ways to check the validity of shower calculations, and 
of the corresponding interaction models, is to compare the results for other detectable 
particles, like for instance muons. Muons have a direct link with neutrinos through 
the decay mechanisms, which are rather well known. Unlike neutrinos, not all muons 
produced in atmospheric showers are available at ground level, due to energy loss and 
decay along their path. Therefore, the data taken at different altitudes by means of 
balloon-born detectors are of particular interest. In this respect, the data obtained by the 
CAPRICE experiment are of particular interest and considered among the most reliable by 
the scientific community. We are interested in particular to the data taken at Lynn Lake 
in 1994 and published in 0. In this paper we describe the results of the comparison of our 
FLUKA calculations with this data set, showing that, under the assumption concerning 
the validity of the input primary spectrum adopted so far, a good agreement with these 
muon data is achieved. 

Another neutrino flux calculation has been recently proposed^], with a model capable 
of obtaining a successful comparison with the CAPRICE data. In this case the normal- 
ization of the resulting neutrino flux is different from the FLUKA one, depending on 
neutrino flavor and energy, but with the general trend of producing a 20 - 25% lower 
normalization in the sub-GeV region. According to the authors, when switching from 
their particle production model to the Bartol one the resulting neutrino fluxes obtained 
with their code (CORT) increase by a factor exceeding 50% in the Sub-GeV range. They 
conclude that their model, which is able to reproduce the muon fluxes in the atmosphere, 
predicts significantly lower neutrino fluxes than those currently used in the experimental 
analyses (mostly the Bartol and IIKKM0 computations). 

2 Monte Carlo setup 

The FLUKA code |l[] is a general purpose Monte Carlo code for the interaction and 
transport of particles. It is built and maintained with the aim of including the best pos- 
sible physical models in terms of completeness and precision. It contains detailed models 
of electromagnetic, hadron-hadron and hadron-nucleus interactions, covering the range 
from the MeV scale to the many-TeV one. Extensive benchmarking against experimental 
data has been produced (see the references in [Q]). In view of applications for cosmic 
ray physics, we have implemented in the FLUKA environment a 3-Dimensional spheri- 
cal representation of the whole Earth and of the surrounding atmosphere. The latter is 
described by a by a proper mixture of N, O and Ar, arranged in 100 concentric spherical 
shells, whose densities follow the known profile of the USA "standard atmosphere". This 
is the standard choice in neutrino flux calculations, although we are aware that important 
differences exist as a function of latitude and season. For a muon flux calculation the 
specific features of the local geographical site and the actual situation on the days of 
measurement could be relevant for a fully significant comparison with data. We have also 
tried to set the boundaries between the atmospheric shells in such a way to reproduce as 



much as possible the different "ffux weighted average depths" at which the CAPRICE94 
measurements were reported. Considering our discrete structure and the adopted atmo- 
sphere profile, we evaluated a systematic error on the atmospheric depth of less than few 
percent. 

The starting point of any calculation of particle fluxes in the atmosphere is the primary 
flux. In order to allow a direct comparison with previously published computations, we 
used the fit of all-nucleon flux proposed by the Bartol group |^. This flux has also the 
merit that, at least up to about hundred GeV it is quite close to the recent measurements 
of proton primaries performed by BESSf^ and AMS|^. A comparison with these data sets 
is shown in Fig. |l], where the all-nucleon spectrum has been divided in the free-proton 
and bounded-nucleon components. To the bounded-nucleon component, mass groups 



heavier than helium contribute. The sum of such contributions, taken from ref. |]T0 
are represented by a continuous line. Nuclei are treated in the framework of the super- 
position model, namely a nucleus of mass number A, and total energy Eq is considered 
to be equivalent to A nucleons each having energy Eq/A. At present, this is the approach 
followed in all other neutrino calculations, except the case of ref.0, where a Glauber-like 
model is used for the primary nuclear component. 

For the purpose of the present calculation, we have modifled the spectrum in order 
to take into account the solar modulation at the time of the experimental measurement 
(August 8th 1994). In order to do that, we have followed the modeling of ref. [jll|. First, a 



Local Interstellar Spectrum has been derived from the Bartol spectrum at solar minimum, 
and than it has been modulated back to the desired level according to the counting rate 
of the Climax neutron monitor. The effect of solar modulation on all-nucleon spectrum 
is shown in Fig. ^. Primary particles, sampled from the continuous energy spectrum, are 
injected at the top of the atmosphere, at about 100 km of altitude. 

The primary flux is assumed to be uniform and isotropic far from the Earth, and it is 
locally modifled by the effects of the geomagnetic fleld[|l^] and solar modulation. Contrary 
to our previous neutrino flux calculation we have switched on the local magnetic fleld also 
during shower development, in order to allow the bending of charged particle trajectories. 

In the simulation we score the muon flux separately for each muon charge, at each 
boundary crossing between the different atmospheric shells, and at the Earth boundary 
as well. We have tested different acceptance cones to deflne the fluxes, using half-cone 
angles of 8, 12 and 20 degrees. At a few altitudes we have also added a cone of 50 degrees 
of half-aperture. The muon flux is recorded as function of laboratory momentum, in 30 
logarithmic bins between 0.2 and 200 GeV/c. In order to compare our results to the 
experimental data, a binning of calculated fluxes has been performed to reproduce the 
momentum ranges considered in 

We have also performed simulation runs in the 1-Dimensional approximation, that is 
aligning all secondary particles to the primary at each interaction. 



3 Results 

As a first step, we have studied the dependence of simulated fluxes, at different heights, 
as a function of the different aperture cones. As a result we may conclude that up to 20 
degrees, the simulation does not show a significant dependence on the aperture cone at all 
heights. In the 50 degrees cone, instead, a loss of muon flux at low muon momentum can 
be observed with respect to the lower apertures. As an example, in Fig. ^ the negative 
muon flux as a function of momentum for different aperture cones is shown at the nominal 
atmospheric depth of 50.7 g/cm^ (a), corresponding to about 20.6 km of altitude (a) and 
at the nominal depth of 308.9 g/cm^ (b), that is at about 9.08 km of altitude. The loss 
at low muon momentum at large angle is explained by muon energy loss and decay. 

The comparison with the experimental data from CAPRICE is performed using the 
muon flux recorded in the 12 degrees half-aperture cone, the one closer to the apparatus 
acceptance. The results, relative to the full 3-Dimensional simulation runs, are summa- 
rized in Fig. ^ for the negative muon measurements at different depths in atmosphere 
and in Fig. ^ for the positive muons. The numeric values of the simulated fluxes reported 
in these figures are reported in Tables ^ and ^ in the Appendix, for the different detector 
depths. The same data are reported as a function of depth, in different momentum in- 
tervals, in Fig. |] and ^, for negative and positive muons respectively. In order to express 
numerically the level of agreement between data and simulation, we report in table |1| 
the weighted average of ratios of Monte Carlo to experimental data at each atmospheric 
depth, as from Fig. ^ and ^, (thus integrating on all momentum bins). In a similar way, 
table ^ shows the average ratios at each momentum bin (integrating on all atmospheric 
depths). Statistical errors have been used to calculate the weights. 

We consider satisfactory the general level of agreement, also taking into account that 
here we are assuming that no systematic error exist for this set of experimental data. 
Instead, as stated before, a few percent systematics in the calculation should come from 
the atmospheric description. We notice a deficit in the simulation with respect to the data 
at very low momenta at the topmost altitude. In the case of negative muons there could 
be at most a 10% deficit in simulation, in the range 1.23-=-2.0 GeV/c. Here the deficit 
is dominated by low values obtained at the shower depths between 100 and 170 g/cm^, 
although both at 77 and 218 g/cm^ (still close to the shower maximum) a ratio close to 
1.0 is again achieved. For positive muons the maximum deficit is limited to 0.75^0.97 
momentum bin, again dominated by the results in the range between 100 and 170 g/cm^. 
However, we want to stress that, for a given interaction model and calculation scheme, 
these values depend also on the choice of primary all nucleon spectrum. Furthermore we 
also expect that our assumption on the total absence of systematics in the experimental 
data is probably too optimistic. With the present level of statistical error we cannot 
identify systematic problems in the /i'^/;U~ ratio at these altitudes. 

Fig. ^ shows the comparison with the measurement at ground level. Numerical values 
are given in table ^. In this case the flux have been calculated in momentum bins which 
are slightly different with respect to those of experimental data. Therefore, in order 
to evaluate the level of agreement, we calculated the ratios of simulated fluxes to the 



Atmospheric Depth (g/cm^) 


MC/Exp. Data {fi~) 


MC/Exp. Data (/i+) 


3.99 


0.84 ± 0.06 


0.78 ± 0.07 


25.7 


0.97 ± 0.07 


1.00 ± 0.08 


50.7 


1.00 ± 0.05 


1.06 ± 0.07 


77.3 


1.04 ± 0.06 


0.96 ± 0.07 


104.1 


0.94 ± 0.05 


0.91 ± 0.04 


135.3 


0.94 ± 0.05 


0.92 ± 0.06 


173.8 


0.98 ± 0.05 


0.92 ± 0.06 


218.8 


1.03 ± 0.04 


1.11 ± 0.07 


308.9 


0.97 ± 0.06 


0.83 ± 0.07 


463.7 


0.96 ± 0.05 


0.97 ± 0.08 


709 


1.04 ± 0.06 


1.27 ± 0.19 


1000 


0.99 ± 0.03 


0.79 ± 0.04 



Table 1: Weighted average ratios of simulated to experimental data for different atmo- 
spheric depths and integrating on the different momentum ranges. 



Muon momentum (GeV/c) 


MC/Exp. Data (/i") 


MC/Exp. Data (^+) 


0.3-0.53 


1.02 ± 0.05 


0.92 ± 0.04 


0.53-0.75 


0.97 ± 0.03 


0.94 ± 0.04 


0.75-0.97 


1.04 ± 0.05 


0.83 ± 0.03 


0.97-1.23 


1.00 ± 0.05 


0.94 ± 0.05 


1.23-1.55 


0.88 ± 0.04 


1.00 ± 0.06 


1.55-2.0 


0.91 ± 0.04 


0.98 ± 0.06 


2.0-3.2 


0.97 ± 0.04 




3.2-8.0 


1.02 ± 0.03 




8.0-40.0 


0.96 ± 0.05 


— - 



Table 2: Weighted average ratios of simulated to experimental data in different momentum 
bins and integrated on the different atmospheric depths. 



values at the same moinentum (above 0.3 GeV/c) taken from the phenomenological fits 
to experimental data shown as continuous lines in Fig. ^ The distribution of these ratios 
are shown in Fig. ^ Here, the present results could point out a possible deficit in the 
calculated /i"'''s with respect to the fi~ . However, it is again important to remember that 
also the modeling of primary spectrum does affect this result, through the neutron to 
proton ratio which in turn is determined by the relative abundances of nuclei and free 
nucleons. 

In order to understand the possible relevance of the 3-dimensional simulation setup 
we have also analyzed data obtained with the 1-dimensional runs. From this comparison 
we can conclude that differences exist at low muon momenta, where a noticeable excess of 



muons is produced in the 1-Dimensional simulations. As an example, we show in Fig. |TG 



the calculated negative muon flux as a function of atmospheric depth in 4 different mo- 



mentum intervals for the 1-Dimensional and 3-Dimensional simulations. Fig. |TT| shows 
similar results for positive muons. In the lowest momentum bins, the 1-D simulation 
yields a larger muon flux. We attribute the difference between the 1-D and 3-D cases 
to a different path length distribution of muons: the inclusion of kinematic production 
angle, together with particle bending in the geomagnetic field, brings to an increase of 
path length and, as a consequence, to a larger decay probability. This argument has been 
discussed by the Bartol group and other authors in ref. [1^, after that preliminary com- 



parisons of the Bartol 1-D results with CAPRICE data were showing a significant excess 



in the calculated flux in the low momentum region, as also reported in |I5|. However, in 
our opinion it is more probable that the main reason for the excess in that calculation 
was an excess of pion multiplicity. In fact the Bartol group is preparing an improved 
version of their interaction model which seems to yield charged pion multiplicities closer 
to those from the FLUKA modelHTBI. Other differences in Monte Carlo calculations of 



muons in the atmosphere between the 1-D and 3-D case have been already reported in 



ref. |13], where the original interaction model from Bartol was used: however, in this case, 
significant excess of calculated muons as compared to HEAT data was observed at high 
altitude for a wide range of momentum. This reinforces our interpretation on the excess 
of multiplicity in the first Bartol model. 

4 Conclusions 

The results presented here demonstrate that the FLUKA model can reasonably reproduce 
the muon fluxes in the atmosphere as measured by the CAPRICE experiment. Of course, 
the accuracy and reliability of such a comparison depends also on the choice of the atmo- 
spheric model and, mostly, on the assumptions concerning the primary spectrum and the 
geomagnetic cutoff. In particular we are aware that the Bartol all nucleon flux should be 



slightly modified in order to be compatible with the latest data on primaries|jT9|. Further- 
more, one should also consider that also the present data from balloon experiments might 
be affected by systematic errors which could be ^ 1O%[0. In this work we are are relying 
on the choices adopted for our calculation of neutrino fluxes in atmosphere, and therefore 



we can use the results of the present study to at least partially validate the neutrino 
results, due to the strong link between muon and neutrino production. As a result we 
can conclude that, having used the same primary spectrum, the normalization difference 
in neutrino flux with respect to the Bartol results cannot be attributed to an insufficient 
particle production yield in the shower development in the FLUKA model. As far as 
absolute normalization is concerned, the muon data considered here allow to conclude 
that the FLUKA calculation should have at maximum a 10% error in the neutrino energy 
range up to few GeVs. Of course this last statement is strongly related to the choice of the 
primary spectrum, and it is not possible at this stage to disentangle this factor from the 
feature of the FLUKA interaction model, which is otherwise independently benchmarked 
against particle production data. We wish to stress that data taken in controlled condi- 
tions at accelerators should remain the main constraint for any interaction model to be 
used for particle production in atmosphere. In particular we remind that hadron-Nucleus 
and Nucleus-Nucleus interactions in the energy range l-i-30 GeV are still insufficiently 
known to allow the construction of a fully reliable model. The check with measured fluxes 
of different particles in atmosphere remains however relevant in order to verify that all 
other ingredients of the simulation setup are taken into account correctly. 

Another conclusion of this work is that 3-Dimensional effects are of some important 
also for low energy muons. If we assume that our 3-D vs. 1-D effect is correct, there could 
be implications for the results recently presented by Fiorentini et al., in 0. However, the 
significant differences in SubGev neutrino fluxes between the results of and those from 
the FLUKA calculation [§] seem hard to explain on the basis of a geometrical effect or of 
differences in the primary spectrum. The main difference should be the interaction model, 
but then it is still difficult to explain why different neutrino fluxes arise from calculations 
yielding reasonable reproduction of the same data on muon fluxes, a part from technical 
errors or, maybe, from important differences in muon energy loss and decay, but at this 
stage this is just a guess. As far as FLUKA is concerned the energy loss models have 
been obviously widely cross-checked and verifled after many years of practice in preci- 
sion problems like radiation shielding, dosimetry and other radio-biological applications. 
Muon decay with correct polarization effects has been checked together with the Bartol 
group. 
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Figure 1: All- nucleon primary spectrum (from Bartol group) used for the FLUKA calcula- 
tions in the atmosphere as compared to recent AMS and BESS data sets. The continuous 



line representing higher mass components has been derived from ref. |^ 
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Figure 2: Ratio of all-nucleon primary spectrum modulated for the date of August S*'' 
1994 to that at solar minimum value, as obtained using the algorithm of ref. [Ill 
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Figure 3: Calculated Negative muon flux as a function of momentum for different aperture 
cones is shown at the atmospheric depth of 50.7 g/cm^ (a), corresponding to about 20.6 
km of altitude (a) and at the depth of 308.9 g/cm^ (b), that is at about 9.08 km of 
altitude). The lines are phenomenological flts meant to guide the eye. 
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Figure 4: Comparison between simulated and detected negative muon flux as a function 
of momentum for different atmospheric depths. 
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Figure 5: Comparison between simulated (small full symbols) and detected positive muon 
flux (open symbols) as a function of momentum for different atmospheric depths. 
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Figure 8: Comparison between simulated (full symbols) and detected (open symbols) 
negative and positive muon fluxes as a function of momentum at ground level. The 
continuous lines are phenomenological fits to the experimental data (above 0.3 GeV/c). 
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Figure 9: Distribution of the ratios of calculated to experimental muon flux at ground 
level for the different muon charges. 
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Figure 10: Comparison between 3-Dimensional and 1-Dimensional negative muon flux 
calculations as a function of atmospheric depth in a few different momentum ranges. The 
lines are drawn to guide the eye. case. 
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Figure 11: Comparison between 3-Dimensional and 1-Dimensional positive muon flux 
calculations as a function of atmospheric depth in a few different momentum ranges. The 
hues are drawn to guide the eye. 
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Appendix: Tables of numerical values 

Table 3: Simulated negative muon flux as a function of momentum for the different depths 
in atmosphere. 



Muon Mom. 


Flux 


Flux 


Flux 


(GeV/c) 


(cm s sr GeV/c)^l 


(m^ ssr GeV/c)-l 


(m2 ssr GeV/c)-l 




3.9 g/cm^ 


25.7 g/cm2 


50.7 g/cm^ 


0.22 


9.31 -10-^ ±36 % 


2.63 • 10-3 ± 19 % 


7.93 • 10-3 ± 14 % 


0.28 


5.18 


10-^ ± 31 % 


6.04 


10-3 ± 24 % 


6.94 


10-3 ± 15 % 


0.36 


1.92 


10-3 ± 25 % 


4.55 


10-3 ± 24 % 


9.45 


10-3 ± 10 % 


0.45 


8.63 


10-^ ± 25 % 


2.85 


10-3 ± 18 % 


8.75 


10-3 ± 12 % 


0.56 


1.17 


10-3 ± 17 % 


2.14 


10-3 ± 17 % 


4.68 


10-3 ± 12 % 


0.71 


6.33 


10-^ ± 22 % 


1.85 


10-3 ± 12 % 


5.15 


10-3 ± 13 % 


0.89 


8.43 


10-^ ± 35 % 


1.68 


10-3 ± 12 % 


3.23 


10-3 ± 12 % 


1.12 


4.47 


10-^ ± 32 % 


9.58 


10-4 ± 20 % 


2.26 


10-3 ± 13 % 


1.42 


2.52 


10-^ ± 18 % 


6.83 


10-4 ± 18 % 


1.70 


10-3 ± 16 % 


1.78 


1.94 


10-^ ± 24 % 


3.84 


10-4 ± 19 % 


1.27 


■ 10-3 ± 9 % 


2.24 


1.10 


10-^ ± 21 % 


1.34 


10-4 ± 23 % 


6.74 


10-4 ± 12 % 


2.83 


7.03 


10-5 ± 37 % 


8.21 


10-5 ± 21 % 


2.46 


10-4 ± 16 % 


3.56 


2.72 


10-5 ± 24 % 


9.20 


10-5 ± 22 % 


1.84 


10-4 ± 15 % 


4.48 


2.18 


10-5 ± 30 % 


5.22 


10-5 ± 25 % 


1.55 


10-4 ± 15 % 


5.64 


6.62 


10-6 ± 37 % 


2.66 


10-5 ± 18 % 


6.18 


10-5 ± 16 % 


7.10 


5.66 


10-6 ± 26 % 


1.29 


10-5 ± 23 % 


5.04 


10-5 ± 23 % 


8.93 


1.13 


10-5 ± 40 % 


1.52 


10-5 ± 30 % 


3.44 


10-5 ± 15 % 


11.25 


1.35 


10-6 ± 40 % 


3.83 


10-6 ± 29 % 


1.00 


10-5 ± 22 % 


14.16 


3.46 


10-^ ± 64 % 


2.11 


10-6 ± 30 % 


5.93 


10-6 ± 19 % 


17.83 


3.90 


10-^ ± 37 % 


1.34 


10-6 ± 23 % 


2.65 


10-6 ± 17 % 


22.44 


1.37 


10-^ ± 60 % 


6.22 


10-^ ± 26 % 


3.11 


10-6 ± 24 % 


28.25 


5.71 


10-s ± 63 % 


4.84 


10-^ ± 22 % 


1.08 


10-6 ± 30 % 


35.57 


7.24 


10-^ ± 67 % 


8.17 


10-^ ± 59 % 


1.26 


10-6 ± 39 % 


44.77 


1.12 


10-^ ± 54 % 


1.56 


10-^ ± 32 % 


2.41 


10-^ ± 24 % 


56.37 


3.62 


10"^ ± 49 % 


6.07 


10-s ± 46 % 


9.43 


10-8 ± 38 % 


70.96 


5.63 


10-9 ± 52 % 


2.97 


10-s ± 37 % 


1.06 


10-^ ± 31 % 


89.34 


5.20 -10-1° ±65 % 


1.14 


10-9 ± 35 % 


8.23 


10-9 ± 55 % 


112.47 


3.95 -10-1° ±68 % 


2.47 


10-8 ± 66 % 


3.74 


10-8 ± 47 % 


141.59 


5.13 -10-11 ±68 % 


1.52 


10-8 ±61 % 


1.56 


10-8 ± 60 % 


178.25 


2.60 -10-1° ±61 % 


1.15 


10-9 ± 69 % 


5.43 


10-9 ± 41 % 




77.3 g/cm^ 


104.1 g/cm^ 


135.3 g/cm^ 


0.22 


9.24 • 10-^ ± 12 % 


1.46 -10"^ ±7% 


1.79 • 10-^ ± 12 % 
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0.28 


1.44 -10-2 ±10 % 


1.64 -10-2 ±10 % 


1.72 • 10-2 ± 10 % 


0.36 


9.90 • 10-3 ± 10 % 


1.78 -10-2 ±9 % 


2.14 -10-2 ±6 % 


0.45 


1.16 -10-2 ±9 % 


1.71 • 10-2 ± 10 % 


1.56 


• 10-2 ± 9 % 


0.56 


9.55 


■ 10-3 ± 8 % 


1.29 • 10-2 ± 8 % 


1.23 


• 10-2 ± 6 % 


0.71 


7.40 


• 10-3 ± 9 % 


1.06 • 10-2 ± 8 % 


9.94 


• 10-3 ± 7 % 


0.89 


6.24 


• 10-3 ± 9 % 


8.84 • 10-3 ± 8 % 


8.47 


• 10-3 ± 7 % 


1.12 


3.67 


• 10-3 ±7% 


6.51 • 10-3 ± 7 % 


6.80 


• 10-3 ± 6 % 


1.42 


2.80 


■ 10-3 ± 8 % 


3.76 • 10-3 ± 10 % 


3.83 


• 10-3 ± 9 % 


1.78 


1.81 


• 10-3 ± 9 % 


2.54 


10-3 ± 10 % 


2.70 


• 10-3 ± 7 % 


2.24 


9.60 • 10-^ ± 11 % 


1.23 


10-3 ± 10 % 


1.50 


• 10-3 ± 8 % 


2.83 


5.56- 


10-4 ± 15 % 


9.43 


10-4 ± 12 % 


1.18- 


10-3 ± 11 % 


3.56 


2.99- 


10-4 ± 14 % 


4.69 


10-4 ± 10 % 


5.32 


• 10-4 ± 8 % 


4.48 


3.35- 


10-4 ± 12 % 


4.66 


10-4 ± 10 % 


5.35 


• 10-4 ± 8 % 


5.64 


1.16- 


10-4 ± 14 % 


1.64 


10-4 ± 10 % 


2.27 


• 10-4 ± 8 % 


7.10 


6.80- 


10-5 ± 18 % 


1.12 


10-4 ± 14 % 


1.39 -10-4 ±11 % 


8.93 


5.13- 


10-5 ± 16 % 


6.49 


• 10-5 ± 9 % 


7.27 • 10-5 ± 8 % 


11.25 


2.55- 


10-5 ± 20 % 


3.41 


10-5 ± 17 % 


3.66 • 10-5 ± 16 % 


14.16 


1.06- 


10-5 ± 15 % 


1.53 


10-5 ± 14 % 


1.79- 


10-5 ± 13 % 


17.83 


4.52- 


10-6 ± 18 % 


5.92 


10-6 ± 16 % 


7.16- 


10-6 ± 14 % 


22.44 


3.39- 


10-6 ± 18 % 


3.86 


10-6 ± 17 % 


4.81- 


10-6 ± 17 % 


28.25 


2.29- 


10-6 ± 21 % 


3.49 


10-6 ± 20 % 


3.86- 


10-6 ± 19 % 


35.57 


2.06- 


10-6 ± 25 % 


2.42 


10-6 ± 21 % 


2.68- 


10-6 ± 21 % 


44.77 


3.88- 


10-^ ± 17 % 


5.05 


10-^ ± 17 % 


5.54- 


10-^ ± 16 % 


56.37 


2.25- 


10-^ ± 31 % 


3.93 


10-^ ± 25 % 


4.71- 


10-^ ± 22 % 


70.96 


1.96- 


10-^ ± 24 % 


2.39 


10-^ ± 21 % 


2.43- 


10-^ ± 21 % 


89.34 


3.37- 


10-s ± 29 % 


4.87 


10-^ ± 24 % 


6.65- 


10-^ ± 22 % 


112.47 


3.98- 


10-8 ± 44 % 


4.48 


10-s ± 39 % 


6.10- 


10-s ± 32 % 


141.59 


2.02- 


10-s ± 51 % 


3.00 


10-s ± 36 % 


3.53- 


10-s ± 32 % 


178.25 


6.61- 


10-9 ± 34 % 


7.53 


10-9 ± 32 % 


9.32- 


10-9 ± 26 % 




173.8 g/cm^ 


218.8 g/cm2 


308.9 g/cm2 


0.22 


1.67 -10-2 ±12 % 


1.53 • 10-2 ± 11 % 


1.88 • 10-2 ± 12 % 


0.28 


1.80 -10-2 ±8 % 


2.16 -10-2 ±8 % 


2.12 -10-2 ±6 % 


0.36 


1.93 


• 10-2 ± 8 % 


2.16 


• 10-2 ± 6 % 


2.10 


• 10-2 ± 9 % 


0.45 


1.85 


• 10-2 ±7% 


1.82 


• 10-2 ± 8 % 


1.91 


• 10-2 ± 8 % 


0.56 


1.41 


■ 10-2 ± 9 % 


1.64 


• 10-2 ± 7 % 


1.48 


• 10-2 ± 6 % 


0.71 


1.22 


• 10-2 ± 6 % 


1.34 


• 10-2 ± 7 % 


1.38 


• 10-2 ± 6 % 


0.89 


9.97 


• 10-3 ± 6 % 


1.03 


• 10-2 ± 7 % 


9.03 


• 10-3 ± 7 % 


1.12 


7.34 


• 10-3 ± 6 % 


6.54 


• 10-3 ± 7 % 


6.77 


• 10-3 ± 8 % 


1.42 


4.26 


• 10-3 ± 9 % 


4.94 


• 10-3 ± 7 % 


4.68 


• 10-3 ± 5 % 
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1.78 


2.77 • 10-3 ± 8 % 


2.59 • 10-3 ± 9 % 


3.07 • 10-3 ± 8 % 


2.24 


1.66-10-3±7% 


1.86 • 10-3 ± 8 % 


1.86 -10-3 ±6 % 


2.83 


1.27 -10-3 ±10 % 


1.27 • 10-3 ± 10 % 


1.36-10-3±7% 


3.56 


6.36 • 10-^ ± 8 % 


7.05 • 10-4 ± 9 % 


8.80 • 10-4 ± 8 % 


4.48 


6.03 • 10-^ ± 7 % 


5.28 • 10-4 ± 6 % 


5.60 • 10-4 ± 6 % 


5.64 


2.58 • 10-^ ± 7 % 


3.07 • 10-4 ± 8 % 


3.21 • 10-4 ± 7 % 


7.10 


1.57 -10-4 ±11 % 


1.67 • 10-4 ± 10 % 


1.99 -10-4 ±10 % 


8.93 


8.34 • 10-5 ± 8 % 


8.62 • 10-5 ± 8 % 


9.93 • 10-5 ± 7 % 


11.25 


3.91 • 10-5 ± 15 % 


4.72 • 10-5 ± 12 % 


6.14 -10-5 ±10 % 


14.16 


1.99 • 10-5 ± 12 % 


2.56 • 10-5 ± 10 % 


2.56 • 10-5 ± 11 % 


17.83 


1.07 -10-5 ±14 % 


1.00 • 10-5 ± 13 % 


1.24 -10-5 ±13 % 


22.44 


5.32 • 10-6 ± 15 % 


5.95 • 10-6 ± 14 % 


7.37 • 10-6 ± 13 % 


28.25 


4.09 • 10-6 ± 18 % 


5.50 • 10-6 ± 14 % 


5.94 • 10-6 ± 13 % 


35.57 


3.45 • 10-6 ± 21 % 


3.58 • 10-6 ± 20 % 


3.98 • 10-6 ± 18 % 


44.77 


6.78 • 10-7 ± 15 % 


7.14 • 10-7 ± 14 % 


9.08 • 10-7 ± 14 % 


56.37 


6.09 • 10-7 ± 19 % 


6.39 • 10-7 ± 19 % 


7.44 • 10-7 ± 17 % 


70.96 


2.78 • 10-7 ± 20 % 


3.10 • 10-7 ± 18 % 


3.46 • 10-7 ± 19 % 


89.34 


9.49 • 10-s ± 28 % 


9.89 • 10-s ± 27 % 


1.08 -10-7 ±27% 


112.47 


6.57 • 10-^ ± 29 % 


7.03 • 10-^ ± 27 % 


7.77 • 10-^ ± 27 % 


141.59 


3.61 -10-^ ±31 % 


4.09 • 10-^ ± 28 % 


4.62 • 10-^ ± 27 % 


178.25 


1.10 -10-^ ±23 % 


1.31 -10-^ ±23 % 


1.45 -10-^ ±22 % 




463.7 g/cm^ 


709.0 g/cm^ 


1000.0 g/cm^ 


0.22 


7.10 -10-3 ±15 % 


9.52 • 10-3 ± 21 % 


2.32 • 10-3 ± 13 % 


0.28 


1.09 • 10-2 ± 14 % 


8.44 • 10-3 ± 9 % 


4.00 • 10-3 ± 23 % 


0.36 


7.59 • lO--'^ ± 10 % 


9.48 • 10-3 ± 10 % 


4.40 • 10-3 ± 16 % 


0.45 


9.05 • 10-3 ± 9 % 


6.95 • 10-3 ± 10 % 


3.54 • 10-3 ± 14 % 


0.56 


8.33 


10-3 ± 6 % 


7.92 • 10-3 ± 7 % 


2.56 • 10-3 ± 13 % 


0.71 


5.87 


10-3 ± 8 % 


6.06 • 10-3 ± 9 % 


2.39 • 10-3 ± 12 % 


0.89 


5.25 


10-3 ± 8 % 


4.59 • 10-3 ± 6 % 


1.91 • 10-3 ± 11 % 


1.12 


3.56 


10-3 ± 7 % 


3.62 • 10-3 ± 7 % 


1.34-10-3±7% 


1.42 


2.40 


10-3 ± 9 % 


2.16 • 10-3 ± 10 % 


1.98 • 10-3 ± 13 % 


1.78 


2.16 


10-3 ± 9 % 


2.11 -10-3 ±8 % 


1.26-10-3±7% 


2.24 


1.80 


10-3 ± 8 % 


1.69 • 10-3 ± 8 % 


9.97 • 10-4 ± 9 % 


2.83 


9.90 


10-^ ± 7 % 


1.03 • 10-3 ± 8 % 


6.98 • 10-4 ± 8 % 


3.56 


7.67 


10-^ ± 8 % 


6.86 • 10-4 ± 7 % 


4.77 • 10-4 ± 10 % 


4.48 


3.62 


10-^ ± 8 % 


3.57 • 10-4 ± 8 % 


3.24 • 10-4 ± 9 % 


5.64 


2.99 


10-^ ± 5 % 


3.00 • 10-4 ± 5 % 


2.13 -10-4 ±10 % 


7.10 


1.39 


10-^ ± 8 % 


1.31 • 10-4 ± 9 % 


1.10 -10-4 ±9 % 


8.93 


8.61 


10-5 ± 8 % 


8.50 • 10-5 ± 8 % 


6.70 • 10-5 ± 9 % 
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11.25 


4.88 • 10-5 ± 9 % 


4.81 • 10-5 ± 9 % 


4.23 • 10-5 ± 12 % 


14.16 


2.72 • 10-5 ± 10 % 


2.71 • 10-5 ± 10 % 


2.57 


10-5 ± 11 % 


17.83 


1.44 -10-5 ±12 % 


1.41 • 10-5 ± 12 % 


1.31 


10-5 ± 13 % 


22.44 


1.08 -10-5 ±9 % 


1.08 • 10-5 ± 9 % 


1.05 


10-5 ± 12 % 


28.25 


7.34 • 10-6 ± 13 % 


7.37 • 10-6 ± 13 % 


6.53 


10-6 ± 12 % 


35.57 


3.50 • 10-6 ± 13 % 


3.47 • 10-6 ± 13 % 


3.41 


10-6 ± 13 % 


44.77 


1.29 • 10-6 ± 13 % 


1.40 • 10-6 ± 14 % 


1.35 


10-6 ± 14 % 


56.37 


9.40 • 10-^ ± 15 % 


8.98 • 10-^ ± 16 % 


9.15 


10-^ ± 15 % 


70.96 


3.27 • 10-^ ± 18 % 


3.27 • 10-^ ± 18 % 


3.32 


10-^ ± 18 % 


89.34 


1.67 -10-^ ±23 % 


1.71 -10-^ ±23 % 


1.66 


10-^ ± 24 % 


112.47 


1.02 -10-^ ±21 % 


1.03 -10-^ ±21 % 


1.00 


10-^ ± 22 % 


141.59 


5.60 • 10-^ ± 22 % 


5.61 • 10-^ ± 22 % 


5.15 


10-^ ± 23 % 


178.25 


1.74 -10-^ ±19 % 


1.76 • 10-^ ± 19 % 


1.75 


10-^ ± 19 % 
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Table 4: Simulated positive muon flux as a function of momentum for the different depths 
in atmosphere. 



Muon Mom. 


Flux 


Flux 


Flux 


(GeV/c) 


(cm2 ssr GeV/c)~l 


(m^ ssr GeV/c)-l 


(m2 ssr GeV/c)-l 




3.9 g/cm2 


25.7 g/cm2 


50.7 g/cm2 


0.22 


8.18 -10-^ ±50 % 


7.37 ■ 10-3 ± 18 % 


1.47 • 10-2 ± 9 % 


0.28 


7.41 


10-^ ± 27 % 


1.08 • 10-2 ± 14 % 


1.09 ■ 10-2 ± 13 % 


0.36 


9.79 


10-^ ± 27 % 


9.51 ■ 10-3 ± 11 % 


1.34 -10-2 ±10 % 


0.45 


1.98 


10-3 ± 32 % 


9.67 -10-3 ±12 % 


1.58 -10-2 ±7% 


0.56 


6.18 


10-^ ± 21 % 


7.39 • 10-3 ± 7 % 


1.31 


• 10-2 ± 7 % 


0.71 


2.05 


10-3 ± 19 % 


7.52 ■ 10-3 ± 12 % 


1.15 


• 10-2 ± 7 % 


0.89 


8.24 


10-^ ± 18 % 


4.18 -10-3 ±10 % 


7.53 


• 10-3 ± 7 % 


1.12 


9.47 


10-^ ± 21 % 


2.81 • 10-3 ± 8 % 


5.07 


• 10-3 ± 8 % 


1.42 


6.02 


10-^ ± 20 % 


2.11 -10-3 ±8 % 


4.70 


• 10-3 ± 9 % 


1.78 


2.37 


10-^ ± 28 % 


1.81 -10-3 ±10 % 


2.18 -10-3 ±10 % 


2.24 


2.20 


10-^ ± 32 % 


8.33 


10-^ ± 16 % 


1.51 


10-3 ± 14 % 


2.83 


7.05 


10-5 ± 23 % 


5.32 


10-4 ± 14 % 


7.43 


10-4 ± 10 % 


3.56 


8.84 


10-5 ± 42 % 


3.88 


10-^ ± 12 % 


5.60 


10-4 ± 12 % 


4.48 


4.09 


10-5 ± 35 % 


1.50 


10-^ ± 13 % 


2.88 


10-4 ± 12 % 


5.64 


8.99 


10-6 ± 39 % 


8.28 


10-5 ±11 % 


1.43 


10-4 ± 10 % 


7.10 


1.18 


10-5 ± 50 % 


4.28 


10-5 ± 15 % 


8.14 


10-5 ± 13 % 


8.93 


3.02 


10-6 ± 35 % 


1.82 


10-5 ±21 % 


4.00 


10-5 ± 18 % 


11.25 


2.57 


10-6 ± 33 % 


1.35 


10-5 ± 19 % 


1.85 


10-5 ± 17 % 


14.16 


1.34 


10-6 ± 32 % 


4.31 


10-6 ± 22 % 


7.12 


10-6 ± 18 % 


17.83 


5.71 


10-s ± 37 % 


6.60 


10-6 ±31 % 


1.05 


10-5 ± 20 % 


22.44 


5.31 


10-^ ± 29 % 


3.30 


10-6 ± 18 % 


4.23 


10-6 ± 15 % 


28.25 


8.30 


10-^ ± 53 % 


1.32 


10-6 ± 18 % 


1.88 


10-6 ± 17 % 


35.57 


8.85 


10-9 ± 37 % 


9.68 


10-^ ± 25 % 


1.44 


10-6 ± 19 % 


44.77 


1.39 


10-s ± 46 % 


3.09 


10-^ ± 30 % 


5.67 


10-^ ± 22 % 


56.37 


2.90 


10-s±51 % 


1.14 


10-^ ± 36 % 


3.83 


10-7 ^ 44 % 


70.96 


4.52 


10-s ± 57 % 


7.48 


10-8 ± 36 % 


1.21 


10-^ ± 26 % 


89.34 


1.09 


10-8 ± 84 % 


1.07 


10-^ ± 64 % 


1.33 


10-^ ± 51 % 


112.47 


1.77 


10-9 ± 77 % 


3.24 


10-8 ± 35 % 


5.60 


10-8 ± 28 % 


141.59 


0.00 • 10° ± % 


3.21 


10-9 ± 42 % 


6.85 


10-9 ± 28 % 


178.25 


0.00 • 10° ± % 


1.99 


10-9 ± 59 % 


2.43 


10-9 ± 48 % 




77.3 g/cm^ 


104.1 g/cm^ 


135.3 g/cm2 


0.22 


1.76 -10-2 ±11 % 


1.91 -10-2 ±9 % 


1.99 -10-2 ±12 % 


0.28 


1.99 -10-2 ±9 % 


2.23 • 10-2 ± 8 % 


2.46 • 10-2 ± 8 % 


0.36 


2.20 ■ 10-2 ± 6 % 


2.47 • 10-2 ± 8 % 


2.51 -10-2 ±8% 
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0.45 


1.84-10-2±7% 


2.49 • 10-2 ± 5 % 


2.13 • 10-2 ± 5 % 


0.56 


1.71 


■ 10-2 ± 9 % 


1.76 


• 10-2 ± 6 % 


2.24 


• 10-2 ± 7 % 


0.71 


1.41 


• 10-2 ± 6 % 


1.86 


• 10-2 ± 6 % 


1.76 


• 10-2 ± 7 % 


0.89 


1.19 


• 10-2 ± 7 % 


1.07 


• 10-2 ± 6 % 


1.22 


• 10-2 ± 6 % 


1.12 


8.09- 


10-3 ± 10 % 


9.85 


• 10-3 ± 6 % 


1.06 


• 10-2 ± 7 % 


1.42 


5.78 


• 10-3 ± 6 % 


7.45 


• 10-3 ± 8 % 


6.54 


• 10-3 ± 7 % 


1.78 


3.34 


• 10-3 ± 9 % 


4.26 


• 10-3 ± 9 % 


4.70 


• 10-3 ± 8 % 


2.24 


2.17 


• 10-3 ± 9 % 


2.62 


• 10-3 ± 8 % 


2.80 • 10-3 ± 10 % 


2.83 


1.11 


• 10-3 ± 7 % 


1.32 


• 10-3 ± 8 % 


1.24 -10-3 ±10 % 


3.56 


6.74 


• 10-^ ± 8 % 


9.12 


• 10-4 ± 7 % 


1.03-10-3±7% 


4.48 


3.92 • 10-^ ± 10 % 


4.13 -10-^ ±10 % 


4.70 • 10-4 ± 9 % 


5.64 


1.95 -10-^ ±9 % 


2.15 -10-4 ±9 % 


2.67 • 10-4 ± 8 % 


7.10 


1.30-10~^±14 % 


1.61 • 10-^ ± 12 % 


1.88 -10-4 ±10 % 


8.93 


5.77- 


10-5 ± 14 % 


6.88- 


10-5 ± 12 % 


8.23- 


10-5 ± 11 % 


11.25 


2.46- 


10-5 ± 15 % 


2.63- 


10-5 ± 12 % 


3.10- 


10-5 ± 11 % 


14.16 


2.65- 


10-5 ± 21 % 


3.35- 


10-5 ± 20 % 


3.77- 


10-5 ± 16 % 


17.83 


1.33- 


10-5 ± 16 % 


1.35- 


10-5 ± 15 % 


1.67- 


10-5 ± 14 % 


22.44 


6.27- 


10-6 ± 19 % 


9.10- 


10-6 ± 14 % 


1.03- 


10-5 ± 13 % 


28.25 


3.59- 


10-6 ± 16 % 


3.94- 


10-6 ± 15 % 


4.36- 


10-6 ± 15 % 


35.57 


1.81- 


10-6 ± 19 % 


3.24- 


10-6 ± 21 % 


3.27- 


10-6 ± 21 % 


44.77 


8.21- 


10-^ ± 21 % 


9.28- 


10-^ ± 20 % 


9.98- 


10-^ ± 18 % 


56.37 


4.88- 


10-^ ± 35 % 


9.58- 


10-^ ± 34 % 


1.05- 


10-6 ± 30 % 


70.96 


1.43- 


10-^ ± 22 % 


1.91- 


10-^ ± 21 % 


2.19- 


10-^ ± 19 % 


89.34 


1.59- 


10-^ ± 43 % 


1.74- 


10-^ ± 39 % 


1.92- 


10-^ ± 35 % 


112.47 


6.66- 


10-^ ± 24 % 


6.99- 


10-^ ± 22 % 


7.90- 


10-^ ± 20 % 


141.59 


1.28- 


10-s ± 21 % 


1.38- 


10-s ± 21 % 


1.74- 


10-8 ± 17 % 


178.25 


3.13- 


10-9 ± 38 % 


3.63- 


10-9 ± 33 % 


3.93- 


10-9 ± 30 % 




173.8 g/cm^ 


218.8 g/cm2 


308.9 g/cm2 


0.22 


2.41 • 10-2 ± 11 % 


2.01 


• 10-2 ± 9 % 


1.63-10-2±7% 


0.28 


2.32 • 10-2 ± 8 % 


2.28 


• 10-2 ± 7 % 


1.56 • 10-2 ± 10 % 


0.36 


2.76 


• 10-2 ± 5 % 


2.03- 


10-2 ± 10 % 


1.55-10-2±7% 


0.45 


2.12 


• 10-2 ±7% 


2.38 


• 10-2 ± 6 % 


2.03 


• 10-2 ± 8 % 


0.56 


2.30 


■ 10-2 ± 8 % 


1.95 


• 10-2 ± 5 % 


1.38 


• 10-2 ± 4 % 


0.71 


1.65 


• 10-2 ± 5 % 


1.46 


• 10-2 ± 6 % 


1.35 


• 10-2 ± 6 % 


0.89 


1.07 


• 10-2 ± 6 % 


1.45 


• 10-2 ± 7 % 


1.09 


• 10-2 ± 5 % 


1.12 


1.14 


• 10-2 ± 5 % 


9.80 


• 10-3 ± 5 % 


7.66 


• 10-3 ± 7 % 


1.42 


5.82 


• 10-3 ± 8 % 


6.22 


• 10-3 ± 5 % 


5.90 


• 10-3 ± 7 % 


1.78 


4.80 


• 10-3 ± 7 % 


4.81 


• 10-3 ± 6 % 


4.44 


• 10-3 ± 9 % 


2.24 


2.68 


• 10-3 ± 7 % 


2.36 


• 10-3 ± 5 % 


2.08 


• 10-3 ± 9 % 
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2.83 


1.47 -10-3 ±9 % 


1.72 • 10-3 ± 7 % 


1.43-10-3±7% 


3.56 


1.10 -10-3 ±6 % 


9.81 • 10-4 ± 8 % 


9.55 • 10-4 ± 7 % 


4.48 


5.58 • 10-^ ± 10 % 


5.23 • 10-4 ± 9 % 


5.03 • 10-4 ± 7 % 


5.64 


3.08 • 10-^ ± 7 % 


3.30 • 10-4 ± 7 % 


3.23 • 10-4 ± 9 % 


7.10 


1.91 -10-^ ±9 % 


2.09 • 10-4 ± 7 % 


2.31 • 10-4 ± 8 % 


8.93 


8.20 • 10-5 ± 12 % 


9.27 • 10-5 ± 10 % 


1.07 -10-4 ±8 % 


11.25 


4.10 • 10-5 ± 10 % 


4.26 • 10-5 ± 10 % 


5.97 • 10-5 ± 7 % 


14.16 


4.17 -10-5 ±14 % 


4.78 • 10-5 ± 13 % 


4.84 • 10-5 ± 10 % 


17.83 


2.35 • 10-5 ± 11 % 


2.61 • 10-5 ± 11 % 


2.97 • 10-5 ± 11 % 


22.44 


1.11 -10-5 ±12 % 


1.22 • 10-5 ± 11 % 


1.21 • 10-5 ± 10 % 


28.25 


4.86 • 10-6 ± 14 % 


5.48 • 10-6 ± 14 % 


5.99 • 10-6 ± 12 % 


35.57 


4.24 • 10-'^ ± 18 % 


5.17 -10-6 ±16 % 


5.36 • 10-6 ± 16 % 


44.77 


1.16 -10-6 ±17% 


1.20 • 10-6 ± 16 % 


1.28 • 10-6 ± 15 % 


56.37 


1.07 -10-6 ±30 % 


1.21 -10-6 ±26 % 


1.29 -10-6 ±25 % 


70.96 


2.49 • 10-^ ± 18 % 


3.87 • 10-^ ± 24 % 


3.99 • 10-^ ± 22 % 


89.34 


2.03 • 10-^ ± 33 % 


1.38 • 10-^ ± 16 % 


1.77 -10-^ ±17% 


112.47 


8.56 • 10-s ± 19 % 


1.01 • 10-^ ± 16 % 


1.03 -10-^ ±13 % 


141.59 


2.02 • 10-s ± 15 % 


2.48 • 10-s ± 16 % 


3.41 • 10-s ± 16 % 


178.25 


4.51 • 10-9 ± 26 % 


5.38 • 10-9 ± 23 % 


6.58 • 10-9 ± 21 % 




463.7 g/cm2 


709.0 g/cm^ 


1000.0 g/cm^ 


0.22 


7.57 • 10-3 ± 12 % 


3.69 • 10-3 ± 26 % 


7.71 • 10-4 ± 23 % 


0.28 


9.69 • 10-3 ± 13 % 


2.27 • 10-3 ± 16 % 


1.47 -10-3 ±28 % 


0.36 


9.22 -10-3 ±11 % 


3.61 • 10-3 ± 11 % 


1.56 • 10-3 ± 18 % 


0.45 


1.01 -10-2 ±9 % 


3.45 • 10-3 ± 17 % 


1.68 -10-3 ±21 % 


0.56 


7.99 • 10-3 ± 8 % 


3.66 • 10-3 ± 11 % 


1.48 • 10-3 ± 16 % 


0.71 


7.65 -10-3 ±11 % 


2.94 • 10-3 ± 8 % 


1.33 • 10-3 ± 16 % 


0.89 


6.60 • 10-3 ± 7 % 


3.09 • 10-3 ± 12 % 


1.02 • 10-3 ± 11 % 


1.12 


4.67 • 10-3 ± 6 % 


2.31 • 10-3 ± 8 % 


1.12 • 10-3 ± 14 % 


1.42 


4.25 • 10-3 ± 7 % 


1.68 -10-3 ±9 % 


1.23 • 10-3 ± 15 % 


1.78 


2.24 • 10-3 ± 5 % 


1.28 • 10-3 ± 11 % 


7.07 • 10-4 ± 12 % 


2.24 


1.91 -10-3 ±9 % 


1.20 • 10-3 ± 9 % 


5.80 • 10-4 ± 12 % 


2.83 


1.28 -10-3 ±6 % 


6.84 • 10-4 ± 9 % 


5.15 -10-4 ±13 % 


3.56 


6.77 • 10-^ ±7% 


5.11 -10-4 ±9 % 


3.68 • 10-4 ± 10 % 


4.48 


4.79 • 10-^ ± 9 % 


3.48 • 10-4 ± 10 % 


3.04 • 10-4 ± 9 % 


5.64 


3.22 • 10-^ ± 8 % 


2.57 • 10-4 ± 8 % 


2.05 • 10-4 ± 8 % 


7.10 


1.95-10-^±7% 


1.42 • 10-4 ± 9 % 


8.94 -10-5 ±11 % 


8.93 


8.71 -10-5 ±9 % 


7.82 • 10-5 ± 10 % 


6.78 • 10-5 ± 9 % 


11.25 


6.17-10-5±7% 


5.93 • 10-5 ± 10 % 


5.92 • 10-5 ± 14 % 


14.16 


5.42 • 10-5 ± 12 % 


4.91 • 10-5 ± 12 % 


3.62 • 10-5 ± 10 % 
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17.83 


2.48 • 10-5 ± 11 % 


2.13 • 10-5 ± 14 % 


2.02 • 10-5 ± 17 % 


22.44 


1.21 


10-5 ± 10 % 


1.13 


10-5 ± 11 % 


9.70 


10-6 ± 13 % 


28.25 


6.08 


10-6 ± 11 % 


7.42 


10-6 ± 13 % 


7.21 


10-6 ± 14 % 


35.57 


5.60 


10-6 ± 20 % 


4.27 


10-6 ± 17 % 


3.96 


10-6 ± 15 % 


44.77 


1.14 


10-6 ± 16 % 


1.15 


10-6 ± 16 % 


9.98 


10-^ ± 15 % 


56.37 


1.35 


10-6 ± 24 % 


1.33 


10-6 ± 24 % 


1.30 


10-6 ± 25 % 


70.96 


4.08 


10-^ ± 22 % 


4.05 


10-^ ± 22 % 


4.06 


10-^ ± 22 % 


89.34 


2.19 


10-^ ± 15 % 


2.70 


10-^ ± 15 % 


2.61 


10-^ ± 16 % 


112.47 


8.56 


10-8 ± 14 % 


8.55 


10-8 ± 15 % 


9.34 


10-8 ± 15 % 


141.59 


3.48 


10-8 ± 16 % 


3.59 


10-8 ± 15 % 


2.88 


10-8 ± 14 % 


178.25 


1.55 


10-8 ± 22 % 


1.56 


10-8 ± 23 % 


1.57 


10-8 ± 22 % 
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Table 5: Simulated muon flux as a function of momentum at ground depth. 



Muon Mom. 


fi'^ Flux 


yU Flux 


(GeV/c) 


(cm2 ssr GeV/c)-l 


(cm^ ssr GeV/c)-l 


0.22 


7.71 • 10^4 ± 23 % 


7.79 • 10-4 ± 24 % 


0.28 


1.47 


10-3 ± 28 % 


8.05 


10-4 ± 17 % 


0.36 


1.56 


10-3 ± 18 % 


1.09 


10-3 ± 20 % 


0.45 


1.68 


10-3 ±21 % 


1.64 


10-3 ± 16 % 


0.56 


1.48 


10-3 ± 16 % 


8.01 


10-4 ± 15 % 


0.71 


1.33 


10-3 ± 16 % 


1.38 


10-3 ± 10 % 


0.89 


1.02 


10-3 ±11 % 


1.29 


10-3 ± 16 % 


1.12 


1.12 


10-3 ± 14 % 


1.26 


10-3 ± 14 % 


1.42 


1.23 


10-3 ± 15 % 


7.66 


10-4 ± 12 % 


1.78 


7.07 


10-4 ± 12 % 


7.35 


10-4 ± 12 % 


2.24 


5.80 


10-4 ± 12 % 


6.77 


10-4 ±11 % 


2.83 


5.15 


10-4 ± 13 % 


4.32 


10-4 ± 12 % 


3.56 


3.68 


10-4 ± 10 % 


3.21 


• 10-4 ± 9 % 


4.48 


3.04 • 10-4 ± 9 % 


2.88 


10-4 ± 10 % 


5.64 


2.05 • 10-4 ± 8 % 


1.31 


10-4 ± 11 % 


7.10 


8.94 -10-5 ±11 % 


9.28 • 10-5 ± 8 % 


8.93 


6.78 • 10-5 ± 9 % 


5.49 • 10-5 ± 9 % 


11.25 


5.92 • 10-5 ± 14 % 


3.68 -10-5 ±11 % 


14.16 


3.62 


10-5 ± 10 % 


1.98 


10-5 ± 13 % 


17.83 


2.02 


10-5 ± 17 % 


1.32 


10-5 ± 12 % 


22.44 


9.70 


10-6 ± 13 % 


9.17 


10-6 ± 13 % 


28.25 


7.21 


10-6 ± 14 % 


5.39 


10-6 ± 15 % 


35.57 


3.96 


10-6 ± 15 % 


3.39 


10-6 ± 14 % 


44.77 


9.98 


10-^ ± 15 % 


1.27 


10-6 ± 14 % 


56.37 


1.30 


10-6 ± 25 % 


8.44 


10-^ ± 15 % 


70.96 


4.06 


10-^ ± 22 % 


3.27 


10"^ ± 18 % 


89.34 


2.61 


10-^ ± 16 % 


1.69 


10-^ ± 24 % 


112.47 


9.34 


10-8 ± 15 % 


9.80 


10-8 ± 21 % 


141.59 


2.88 


10-8 ± 14 % 


5.12 


10-8 ± 23 % 


178.25 


1.57 


10-8 ± 22 % 


1.49 


10-8 ± 19 % 
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